Conformational Transitions of Non-Grafted Polymers Near an Adsorbing Substrate 
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We have performed multicanonical chain-growth simulations of a polymer interacting with an 
adsorbing surface. The polymer, which is not explicitly anchored at the surface, experiences a 
hierarchy of phase transitions between conformations binding and non-binding with the substrate. 
We discuss the phase diagram in the temperature-solubility plane and highlight the transition 
"path" through the free-energy landscape. 

PACS numbers: 05.10.-a, 87.15. Aa, 87.15.Cc 



The recent developments in single molecule experi- 
ments at the nanometer scale, e.g., by means of atomic 
force microscopy (AFM) pj and optical tweezers Q, al- 
low now for a more detailed exploration of structural 
properties of polymers in the vicinity of adsorbing sub- 
strates. The possibility to perform such studies is of es- 
sential biological and technological significance. From 
the biological point of view the understanding of the 
binding and docking mechanisms of proteins at cell mem- 
branes is important for the reconstruction of biological 
cell processes. Similarly, specificity of peptides and bind- 
ing affinity to selected substrates could be of great impor- 
tance for future electronic nanoscale circuits and pattern 
recognition devices. Since single-molecule experiments 
allow monitoring of polymer chains adsorbed at surfaces, 
the investigation of structural deformations of the poly- 
mer shape near substrates is a central aspect of experi- 
mental, computational, and theoretical studies 

In computer simulations and analytical approaches, 
typically, one end of the polymer is anchored at a flat 
substrate and the influence of adhesion and steric hin- 
drance 0, H II H 0, pulling forces [H El or ex- 
ternal fields [12] on the shape of the polymer is con- 
sidered. The question how a flexible substrate, e.g., a 
cell membrane, bends as a reaction of a grafted poly- 
mer, was, for example, addressed in Ref. |l3j . Proteins 
exhibit a strong specificity as the affinity of peptides to 
adsorb at surfaces depends on the amino acid sequence, 
solvent properties, and substrate shape. This was ex- 
perimentally and numerically studied, e.g., for peptide- 
metal [3 Qj5| and peptide-semiconductor [TfjL lrn| inter- 
faces. Binding/folding and docking properties of lattice 
heteropolymers at an adsorbing surface were subject of 
a recent numerical study [l8| . 

In this work we investigate in detail the temperature 
and solubility dependence of adsorption properties for 
a polymer which is not fixed at the surface of the sub- 
strate with one of its ends. This model was inspired by 
the experimental setup used in Refs. [Til where the 
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peptides are initially freely moving in solution before ad- 
sorption. Therefore, there are two main differences in 
comparison with studies of polymers explicitly grafted 
at the substrate: First, the chain can completely desorb 
from the substrate allowing for the investigation of the 
binding/unbinding transition. Second, adsorbed confor- 
mations are possible, where none of the two polymer ends 
is in contact with the surface. 

We use a lattice model 5] for a polymer near an adsorb- 
ing substrate where nearest-neighbor contacts between 
monomers (nonadjacent along the chain) and contacts 
between monomers and the substrate are assigned differ- 
ent energy scales. The polymer energy is given by 

E = -e s n s - s m n m , (I) 

where n s ^ m are the numbers of contacts with the surface 
and between the monomers, respectively. The associated 
energy scales are set to e s — 1 and e m = s, for con- 
venience. The solvent parameter s takes account of the 
goodness of the implicit solvent and rates the two energy 
scales. As a first guess, s > for a poor solvent, where 
globular conformations of the polymer are preferred. For 
s = 1, adsorption and collapse are equally attractive for 
the polymer. Negative values of s refer to good solu- 
bility and stretched conformations dominate. The par- 
tition sum per surface area A for the polymer with the 
described properties at inverse temperature (3 = l/fc^T 
(ks = f in the following) can be written as 

Z((3, s )/A= ]T 9ns n m e^ +sn -\ (2) 

where g na ,n m = ^n s o.g" m + (1 - ^n s o)9n 3 n m is tlie con- 
tact density. In this decomposition, stands for the 
density of unbound conformations, whereas <?„ „ is the 
density of surface and intrinsic contacts of all conforma- 
tions bound to the substrate. Since the number of un- 
bound conformations in the half-space accessible to the 
polymer is unrestricted, formally diverges. For reg- 
ularization, we introduce an impenetrable (but neutral, 
i.e., non-adhesive) wall at a sufficiently large distance z w 
from the substrate, in order to keep its influence on the 
unbound polymer small. 

We apply the powerful multicanonical chain-growth al- 
gorithm, originally introduced for the simulation of lat- 
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FIG. 1: (Color online) Pseudo-phase diagram of a polymer 
with 100 monomers obtained from the specific heat CV as a 
function of temperature T and solubility parameter s. The 
white and black drawn lines indicate the ridges of the specific 
heat profile (see text). The dashed line separates the regions 
of poor and bad solvent (s > 0) from good solvent (s < 0). 



tice proteins [19j, in order to simulate the contact den- 
sity g n „n m directly. With this algorithm, which sets up 
on PERM chain-growth [20|. the contact distributions 
are flattened in a recursive way. This enables optimal 
sampling of the conformational space, and all energetic 
quantities such as, e.g., the specific heat are obtained 
by reweighting the density of contacts with respect to 
temperature and solubility. The main advantage is that 
the whole phase diagram can in principle be constructed 
within a single simulation plj . This method is the key 
for unravelling the detailed structure of the phase di- 
agram, in particular, at low temperatures, where most 
importance sampling Monte Carlo algorithms run into 
difficulties. In order to break correlations being inherent 
in the chain-growth process, we averaged over indepen- 
dent simulations (including the determination of the mul- 
ticanonical weights), and a total statistics of more than 
10 9 chains was accumulated in the production runs for a 
homopolymer with 100 monomers. For confirmation, we 
also investigated polymers with up to 200 monomers [22| • 
In Fig. n we show the pseudo-phase diagram of the 
lOOmer near an attractive substrate and the steric wall in 
a distance z w — 200 from the substrate. The color codes 
the height of the specific heat CV as the function of the 
temperature T and the solvent parameter s; the brighter 
the larger the value of Cv ■ The white and black lines em- 
phasise the maxima of the specific heat which shall serve 
as an orientation for the phase boundaries. While the 
white lines indicate pseudo-transitions being specific for 
the lOOmer, lines drawn in black separate regions which 
are expected to be phases in the strict thermodynamic 
sense. The precise locations of transition lines in the 
thermodynamic limit N — > oo will, however, differ from 
the position for the finite-length system under study. We 
distinguish six thermodynamic phases, four for the ad- 



sorbed (AC1, AC2, AE1, AE2) and two for the desorbed 
(DC, DE) polymer .23]. In the adsorbed-collapsed phase 
AC1, all monomers are in contact with the substrate and 
the two-dimensional (single-layer) conformation ( "film" ) 
is very compact. The transition from AC1 to AC2 is 
the layering phase transition from single to double-layer 
conformations. The white transition lines within AC2 in- 
dicate pseudo-transitions to compact conformations with 
more than two layers. These transitions are expected to 
disappear in the thermodynamic limit 0. The transi- 
tion line between AC1 and the adsorbed-expanded phase 
AE1 is the two-dimensional collapse. It separates the 
compact single-layered conformations in AC1 from the 
dissolved, but still basically two-dimensional conforma- 
tions. White lines in AE1 indicate conformational tran- 
sitions to unstructured conformations extending partially 
into the third dimension. The substrate-contacting layer 
is dissolved and although several layers can form, no ex- 
plicit layering transitions are observed in this region. In 
contrast to AE1, the conformations dominating phase 
AE2 possess a very compact surface layer but less com- 
pact upper layers. As in AE1, the formation of higher- 
order layers is not accompanied with noticeable confor- 
mational transitions. The difference between AE1 and 
AE2 becomes more apparent when approaching the un- 
binding transition line to the phases DE and DC, where 
the polymer has completely desorbed from the substrate: 
In the desorbed-expanded phase DE random-coil confor- 
mations dominate, while in the desorbed-collapsed phase 
DC globular conformations are favored. Phases DE and 
DC are separated by the transition line indicating the 
three-dimensional Q collapse. 

The main difference to a polymer which is explicitly an- 
chored at the substrate with one of its ends, is the occur- 
rence of the strong binding/unbinding transition between 
the A and D phases. In the D phases, the polymer can 
move freely within the cavity, restricted only by the pres- 
ence of the two impenetrable walls. This transition influ- 
ences, however, also the conformational behavior of the 
polymer in the phases AE1 and AE2, the latter not being 
present for the anchored polymer. In fact, phases AC2, 
AE2, and DC lie within the DC/SAG (surface-attached 
globule) regime of the anchored polymer 0, whose 
precise phase structure is not yet completely clarified. 
The phases AE1, AC1, and AC 2 approximately coincide 
in the two systems for low temperatures, although the 
system in our study has more entropic freedom since also 
adsorbed conformations are possible, where none of the 
ends is anchored at the substrate. This may have conse- 
quences for the location of transition lines. 

Apart from the thermodynamic transitions in the tra- 
ditional meaning we also see pseudo-transitions being 
specific to the chosen number of monomers (e.g., a reori- 
entation transition from the 5x5x4 cube with 25 surface 
contacts to the rotated cube exhibiting only 20 contacts 
with the substrate, which occurs within phase AC2 at 
T w 0.7). We are convinced that growing experimental 
capabilities will allow the observation of these effects also 
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FIG. 2: Map of free energy minima identified in the space 
of contact numbers n s with the substrate and n m between 
monomers. Also shown are exemplified "paths" through the 
free-energy landscape for different fixed solvent parameters s. 
The labels I to VII refer to the pseudo phases in the case 
s = 1, described in detail in the text and in Table The lines 
are only guides to the eye. 

for short synthetic or naturally occurring polymers (e.g., 
peptides) . Details of the different conformational phases 
will be reported elsewhere [22| . 

Thermodynamically, the conformations dominating a 
certain phase correspond to the minimum of a suitably 
coarse-grained free energy depending on a few charac- 
teristic observables of the system. The complexity of 
the free energy landscape and its dependence on exter- 
nal parameters such as temperature or solvent strongly 
influences the kinetics of phase transitions. For the poly- 
mer near an adsorbing surface we choose the numbers 
of monomer-substrate contacts, n s , and those between 
monomers, n mi as system observables. According to 
Eq. (J2J, the probability for a polymer conformation is 
given by p(n s ,n m ) oc g n „n m exp([n s + sn m ]/T) and the 
contact free energy reads as 

F s>T {n s ,n m ) = -Tin p(n s ,n m ), (3) 



p{n S) n m ) 




FIG. 3: Probability distribution p(n s ,n m ) for the lOOmer in 
solvent with s = 1 at T = 0.49, where the polymer experiences 
the layering transition from single to double layer. 
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TABLE I: "Path" through the landscape of free energy min- 
ima for a lOOmer in solvent with solubility s — 1 with increas- 
ing temperature. Monomers in contact with the adsorbing 
substrate are shaded in light grey. 



phase T n s n m typical conformations 




where the temperature T and the solubility s are fixed 
external parameters. 

In Fig. [3 we have included all minima of the con- 
tact free energy for the parameter set T £ [0,4.0], 
s G [—2.0,4.0]. Given a fixed solubility s, the stability 
of a conformation with minimal free energy is connected 
with the range of temperatures AT over which the asso- 
ciated free energy is actually the global free energy min- 
imum. We have included into Fig. [21 for several fixed 
solubilities the "paths" of free energy minima hit when 
increasing the temperature from T — up to T — 4, that 
is moving from right to left. 

As an example we consider the case s = 1 for the 
whole region of temperatures. In Table [I] we have listed 
the conformational transitions the lOOmer experiences by 
increasing the temperature. We start at T = with the 
ground state which is a lamellar two-dimensional con- 
formation with 100 surface contacts and 81 monomer- 
monomer contacts. This is the maximal compact con- 
formation that is possible in two dimensions. Note that 
with a conformation characterising a phase we mean all 
structures with the same number of surface and intrin- 
sic contacts as the minimum of the free energy. The 
conformation characterising phase I is highly degenerate. 
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There are in fact about 10 14 (including all symmetries 
except translation) different realizations of the ground 
state. The ground-state conformation remains stable un- 
til T k, 0.2, where the structures become less ordered. 
The lamellar structure is dissolved, but all in all they 
are still two-dimensional and very compact. The confor- 
mational changes in the transition from I to II are rather 
local - in contrast to the probably actual phase transition 
from II to III at T w 0.5, where the number of surface 
contacts is drastically reduced to about half the value of 
the ground-state conformation and thus a second layer 
forms. In Fig. [21 this transition appears as a jump from 
the surface state (n s ,n m ) = (100, 77) to (58, 108). As can 
be seen in Fig. the probability distribution p(n s ,n m ) 
exhibits two distinct peaks at this temperature, which is 
interpreted as strong signal for a first-order transition. 
Entering regime IV, i.e., the adsorbed-expanded phase 
AE2, the dissolution of the surface-contacting layer be- 
gins. This process continues after passing the pseudo- 
transition line to section V, where higher-order layers 
form. Respective bottom layer and upper layers still form 
connected parts - in contrast to phase VI (which belongs 
to AE1), where upper layers can break apart and form 
isolated islands. Increasing the temperature further, we 
approach the second strong first-order-like transition line 
and the polymer unbinds from the substrate for tempera- 
tures T > 2.2. The free energy minimum jumps from the 
contact state (46, 65) to (0, 62) discontinuously (see once 
more Fig. [2J. Therefore, the conformations occurring in 
phase VII do not longer prefer surface contact. Since the 
thermal energy superimposes the relatively weak attrac- 
tion between the monomers at these temperatures, the 
lOOmer in solvent with s = 1 does not experience the 
three-dimensional O transition, because it is already in 



the random coil phase after the unbinding. Note that 
we have also missed the two-dimensional collapse on 
the substrate. For this to happen, the quality of the 
solvent would have to be better (i.e., smaller values of 
s). From the free-energy perspective, both collapse tran- 
sitions are of second order, since the free-energy min- 
ima of the lOOmer at (n a , 100) (two-dimensional surface- 
layer conformations) and (0, n m ) (three-dimensional con- 
formations without contact to the surface) change con- 
tinuously for increasing temperature (see Fig.^J. 

In this work, we have qualitatively analysed the com- 
plete phase diagram for a polymer with 100 monomers 
in solvent near an adsorbing substrate. By means of an 
analysis of the global minima in the free-energy land- 
scape we discussed conformational transitions in the tem- 
perature and solubility parameter space. Two types of 
transitions are experienced by the polymer, phase tran- 
sitions in the thermodynamic sense and transition-type 
cross-over effects which are specific to the given finite 
number of monomers. In the first case, further simula- 
tions of longer polymers combined with finite-size scaling 
analyses will give estimates for the associated transition 
lines. Physically perhaps even more interesting, however, 
are the geometrically induced cross-over effects which are 
expected to become more and more important as the 
high-resolution experimental equipment allows concrete 
measurements in the nanometer range and the design 
of nanoscale devices will take advantage of the specific 
properties of finite- length polymers. 
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